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Description 

[0001] La presente invention a pour objet un procede 
de mesure du temps de propagation d'un signal acous- 
tique dans un ecoulement de fluide entre un premier 
transducteur jouant le role d'emetteur et un second 
transducteur jouant le role de recepteur et situe a une 
distance determinee du premier transducteur, le signal 
acoustique emis par le premier transducteur etant cons- 
titue par au moins une impulsion emise a une frequence 
acoustique determinee Fa et le signal acoustique regu 
par le second transducteur comprenant une serie d'os- 
cillations caracteristiques dont I'amplitude est d'abord 
croissante sur plusieurs periodes, puis decroissante sur 
plusieurs periodes suivantes, I'enveloppe des oscilla- 
tions caracteristiques presentant la forme d'un fuseau, 
le procede consistant a echantillonner le signal acous- 
tique regu, a une frequence d'echantillonnage Fe, a nu- 
meriser le signal acoustique regu echantillonne, et a re- 
chercher, par analyse du signal acoustique regu echan- 
tillonne et numerise, le premier passage a zero signifi- 
catif des oscillations caracteristiques du signal acousti- 
que regu. 

[0002] II est connu depuis de nombreuses annees de 
mesurer le debit d'un fluide (ou son volume) s'ecoulant 
dans une conduite en utilisant la propagation des si- 
gnaux acoustiques emis entre deux transducteurs 
acoustiques situes en des points espaces dans la direc- 
tion de recoupment du fluide. En principe, un signal 
acoustique emis du premier transducteur vers le second 
transducteur est regu par ce second transducteur et le 
temps de propagation Td de ce signal acoustique est 
mesure. De meme, on mesure le temps de propagation 
Tu d'un signal acoustique emis du second transducteur 
vers le premier transducteur apres reception dudit si- 
gnal par ce premier transducteur. 
[0003] Dans un compteur de fluide, le debit peut etre 
obtenu en combinant une mesure des temps de propa- 
gation des deux signaux acoustiques emis entre les 
deux points dans des directions opposees avec une me- 
sure des dephasages acoustiques induits dans chaque 
signal acoustique par la propagation de chacun des si- 
gnaux acoustiques dans recoupment. La demande de 
brevet europeen n° 0 426 309 decrit un exemple d'un 
tel systeme de mesure de debit, dans lequel le signal 
regu est echantillonne et converti numeriquement, la 
mesure du dephasage acoustique etant effectuee en 
realisant une detection synchrone du signal numerise. 
[0004] Lors de la mesure de la vitesse d'ecoulement 
d'un gaz dans un compteur de gaz utilisant deux trans- 
ducteurs ultrasonores, il est necessaire de mesurer le 
temps de parcours de Ponde ultrasonore entre I'instant 
d'emission et I'instant de reception, des lors que la ce- 
lerite des ondes ultrasonores est dependante de la na- 
ture du gaz. 

[0005] On a represents sur la Figure 2 la forme d'un 
signal impulsionnel rectangulaire S 1 de largeur T emis 
a un instant T 0 par un premier transducteur ultrasonore 



dispose dans le flux d'ecoulement d'un fluide en un pre- 
mier point, et la forme du signal S 2 constituant la repon- 
se impulsionnelle regue en echo a un instant T 1 par un 
deuxieme transducteur ultrasonore dispose dans le flux 
5 d'ecoulement du fluide en un deuxieme point distinct du 
premier point. 

[0006] Le signal acoustique S 2 regu par le second 
transducteur est constitue par une serie d'oscillations 
caracteristiques O c qui augmentent d'amplitude sur plu- 

10 sieurs periodes puis decroissent, I'enveloppe des os- 
cillations caracteristiques presentant la forme d'un fu- 
seau. Les oscillations caracteristiques O c du signal S 2 
sont precedees et suivies d'oscillations parasites O p de 
faible amplitude. Pour determiner I'instant T 1 de debut 

15 des oscillations caracteristiques, il convient de reperer 
le premier passage a zero significatif des oscillations ca- 
racteristiques O c du signal acoustique regu S 2 . 
[0007] On a represents sur la Figure 3, a une echelle 
plus grande que sur la Figure 2, un exemple de signal 

20 acoustique S 2 regu en echo a une impulsion rectangu- 
laire S 1 emise a une frequence acoustique determinee 
Fa. 

[0008] Pour determiner le debut des oscillations ca- 
racteristiques O c , selon un procede connu, on fixe une 

25 tension de seuil V s , par rapport a laquelle on compare 
le niveau du signal acoustique regu S 2 , la comparaison 
s'effectuant sur un signal numerique obtenu apres 
echantillonnage du signal acoustique analogique regu, 
a une frequence d'echantillonnage Fe qui constitue par 

so exemple un multiple de la frequence acoustique Fa. 
[0009] Dans ce cas, on repere I'instant T 2 auquel I'am- 
plitude du signal regu franchit la tension de seuil V s , puis 
on identifie I'instant du passage a zero precedent (ou 
suivant) que Ton considere comme I'instant T 1 de depart 

35 des oscillations caracteristiques O c du signal acousti- 
que regu S 2 . 

[0010] Une telle methode de mesure peut conduire a 
des erreurs des lors que les oscillations caracteristiques 
O c du signal acoustique regu S 2 peuvent etre plus ou 

40 moins amplifiees en fonction de la nature du gaz. Ainsi, 
on a represents sur la Figure 4 une courbe S 21 qui cor- 
respond a la forme d'un signal acoustique regu pourde 
I'azote (N 2 ) et une courbe S 22 qui correspond a la forme 
d'un signal acoustique regu pour un melange de gaz car- 

45 bonique et de methane (C0 2 /CH 4 ). On constate que la 
courbe S 21 franchit la tension de seuil V s a un instant 
T 4 qui declenche ('identification du passage a zero pre- 
cedent a un instant T 3 qui est considere a juste titre com- 
me la marque du debut des oscillations caracteristiques 

50 O c . En revanche, on peut constater que la courbe S 22 , 
qui est en phase avec la courbe S 21 , franchit la tension 
de seuil V s a un instant T 6 posterieur a I'instant T 4 et 
decale par rapport a celui-ci de la valeur d'une periode 
T R du signal regu. L'instant T 5 identifie alors pour la 

55 courbe S 22 comme le point de passage a zero precedant 
immediatement le franchissement de seuil a I'instant T 6 , 
est alors considere comme la marque du debut des os- 
cillations caracteristiques O c de la courbe S 22 . Or, com- 
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me on peut le voir sur la Figure 4, la courbe S 2 2 presente 
un lobe negatif qui ne vient qu'affleurer la valeur de la 
tension de seuil V s sans atteindre ou franchir ce seuil. 
[001 1 ] Par suite des attenuations ou amplifications du 
signal regu S 2 , qui varient en fonction de la nature du 
gaz, la methode classique de mesure du temps du pre- 
mier passage a zero des oscillations caracteristiques du 
signal regu peut engendrer une erreur de plus ou moins 
une periode qui fait diminuer de fagon significative la 
precision des mesures. 

[0012] L'invention vise a remedier aux inconvenients 
precites et a permettre de diminuer la sensibilite du pro- 
cede de mesure du temps de passage a zero d'un signal 
acoustique regu, vis-a-vis des variations de formes 
d'onde de ce signal ainsi que vis-a-vis d'eventuelles per- 
turbations exterieures susceptibles d'etre detectees par 
une simple methode de surveillance de franchissement 
de~seuil, et qui pourraient engendrer des mesures erro- 
nees. 

[0013] Ces buts sont atteints, conformement a l'in- 
vention, grace a un procede de mesure du temps de 
propagation d'un signal acoustique dans un ecoulement 
de fluide entre un premier transducteur jouant le role 
d'emetteur et un second transducteur jouant le role de 
recepteur et situe a une distance determinee du premier 
transducteur, le signal acoustique emis par le premier 
transducteur etant constitue par au moins une impulsion 
emise a une frequence acoustique determinee Fa et le 
signal acoustique regu par le second transducteur com- 
prenant une serie d'oscillations caracteristiques dont 
ramplitude est d'abord croissante sur plusieurs perio- 
des, puisdecroissante sur plusieurs periodessuivantes, 
I'enveloppe des oscillations caracteristiques presentant 
la forme d'un fuseau, le procede consistant a echan- 
tillonner le signal acoustique regu, a une frequence 
d'echantillonnage Fe, a numeriser le signal acoustique 
regu echantillonne, et a rechercher, par analyse du si- 
gnal acoustique regu echantillonne et numerise, le pre- 
mier passage a zero significatif des oscillations carac- 
teristiques du signal acoustique regu, caracterise en ce 
que, pour rechercher le premier passage a zero signifi- 
catif des oscillations caracteristiques du signal acousti- 
que regu, on definit au prealable une premiere periode 
caracteristique ideale determinant le premier passage 
par zero des oscillations caracteristiques du signal 
acoustique regu, en caracterisant cette periode ideale 
par un rapport d'amplitudes theorique A entre les ampli- 
tudes maximales Pi- et Pi+ des deux lobes de cette pe- 
riode, on determine pour chaque periode du signal 
acoustique regu echantillonne et numerise, les amplitu- 
des maximales P- et P+ des deux lobes de la periode 
examinee, on compare le rapport de ces amplitudes P- 
et P+ au rapport d'amplitudes theorique correspondant 
de la periode ideale puis, si le resultat de la comparaison 
est superieur a une valeur de seuil G s , on considere la 
periode examinee comme une periode parasite corres- 
pondant a du bruit tandis que, si le resultat de la com- 
paraison est inferieur a cette valeur de seuil G s , on con- 



sidere la periode examinee comme une periode carac- 
teristique et on determine alors le passage par zero en- 
tre deux lobes de cette periode caracteristique, lequel 
passage a zero est considere comme le premier passa- 
5 ge a zero significatif des oscillations caracteristiques du 
signal acoustique regu. 

[0014] Le rapport d'amplitudes theorique A entre les 
amplitudes maximales Pi-et Pi+ des deux lobes de la- 
dite periode ideale est determine au prealable a partir 

10 d'une moyenne, pour plusieurs gaz differents et a diffe- 
rents debits, du rapport entre les amplitudes maximales 
P- et P+ de la periode caracteristique observee a partir 
de I'enregistrement de signaux acoustiques regus. 
[0015] Le rapport d'amplitude A entre les amplitudes 

15 maximales Pi- et P\+ des deux lobes de la premiere pe- 
riode caracteristique ideale d'un signal acoustique regu 
est quasi constant en fonction de la nature des gaz et 
est independant d'un facteur de gain. Le choix d'un cri- 
tere de comparaison base sur ce rapport A permet de 

20 diminuer tres fortement la dependance de la mesure a 
la nature du gaz, et done d'ameliorer la precision et la 
fiabilite des mesures. 

[0016] Selon un premier mode de realisation de l'in- 
vention, la comparaison entre le rapport des amplitudes 
25 P- et P+ de la periode examinee et le rapport d'amplitu- 
des theorique A correspondant de la periode ideale s'ef- 
fectue en calculant un critere de ressemblance G qui 
constitue le resultat de la comparaison compare a la va- 
leur de seuil G s , et qui est defini de la fagon suivante : 

30 



35 



45 



50 



(P-)- A (P+) 



(P-) 



[0017] Selon un deuxieme mode de realisation de in- 
vention, la comparaison entre le rapport des amplitudes 
40 p. et P+ de la periode examinee et le rapport d'amplitu- 
des theorique A correspondant de la periode ideale s'ef- 
fectue en calculant un critere de ressemblance G qui 
constitue le resultat de la comparaison compare a la va- 
leur de seuil G s , et qui est defini de la fagon suivante : 



(P-)- A (P+) 



(P-) -z 



ou z represente un terme de decalage destine a rendre 
le critere de ressemblance G legerement dependant de 
ramplitude maximale (P-) du lobe negatif de la periode 
55 caracteristique examinee. 

[0018] Selon ce deuxieme mode de realisation, on li- 
mite les risques de pouvoir detecter des formes d'ondes 
ressemblant au critere de ressemblance ou gabarit G, 
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qui seraient noyees dans du bruit. 
[0019] Pour reduire encore la probability de detecter 
des periodes ressemblantes de faible amplitude, seion 
un troisieme mode de realisation de Tinvention, la com- 
paraison entre le rapport des amplitudes P- et P+ de la 
periode examinee et le rapport d'amplitudes theorique 
correspondant de la periode ideale s'effectue en calcu- 
lant un critere de ressemblance G qui constitue le resul- 
tat de la comparaison compare a la valeur de seuil G s , 
et qui est defini de la fagon suivante : 



(P-)- A [(P+)-2'] 



ou z represente un terme de decalage destine a rendre 
le critere de ressemblance G legerement dependant de 
I'amplitude maximale (P-) du lobe negatif de la periode 
caracteristique examinee, et z' represente un terme de 
decalage destine a rendre le critere de ressemblance G 
legerement dependant de I'amplitude maximale (P+) du 
lobe positif de la periode caracteristique examinee. 
[0020] Afin d'eviter des erreurs dans la determination 
des amplitudes maximales P- et P+ des deux lobes de 
chaque periode examinee, erreurs qui peuvent etre 
dues a un dephasage existant entre le signal acoustique 
regu et le signal d'echantillonnage, on peut utiliser un 
algorithme d'approximation du maximum (en valeur ab- 
solue) a partir des points echantillonnes au voisinage 
de chaque maximum en vue d'obtenir unemeilleure pre- 
cision sur les valeurs P- et P+. Toujours pour eviter ces 
memes erreurs, lorsque le signal acoustique regu est 
repetable, on effectue plusieurs mesures successives 
respectivement sur plusieurs signaux acoustiques con- 
secutifs, en dephasant pour chaque signal acoustique 
le signal d'echantillonnage par rapport au signal acous- 
tique considere, afin d'obtenir plusieurs points echan- 
tillonnes decales au voisinage du maximum. 
[0021] Selon une variante de realisation permettant 
egalement de reduire la probability de detecter des pe- 
riodes ressemblantes de faible amplitude, on effectue, 
pour chaque periode examinee, deux ou trois calculs du 
critere de ressemblance G avec des valeurs du rapport 
d'amplitudes theorique differentes voisines du rapport 
d'amplitude theorique predetermine et ne s'ecartant pas 
de celle-ci de plus de 10%, et Ton considere la periode 
examinee comme une periode caracteristique si le cri- 
tere de ressemblance calcule G est dans tous les cas 
inferieur a la valeur de seuil G s . 
[0022] Pour determiner la valeur de seuil G s du critere 
de ressemblance G, ainsi que la valeur du terme de de- 
calage z, on trace au prealable, pour differents gaz et 
differents debits, en faisant varier le dephasage entre la 
frequence d'echantillonnage Fe et le signal acoustique 
regu et en appiiquant un facteur de gain variant entre 



0,45 et 1,5, un reseau de premieres courbes represen- 
tant les valeurs maximums du critere de ressemblance 
G pour les periodes caracteristiques observees et un 
reseau de deuxiemes courbes representant les valeurs 
5 minimums du critere de ressemblance G pour les perio- 
des parasites observees, en fonction de differentes va- 
leurs possibles du terme de decalage z, et on choisit la 
valeur de seuil G s et cede du terme de decalage z en 
fonction de la marge de securite possible dans la zone 

io intermediate entre le reseau de premieres courbes et 
le reseau de deuxiemes courbes. 
[0023] A titre d'exemple, la valeur de seuil G s peut 
etre comprise entre 0,7 et 1,7 et la valeur du terme de 
decalage z peut etre comprise entre 0,21 et 0,25 V. 

15 [0024] [.'invention est avantageusement applicable a 
un procede de mesure de la vitesse d'ecoulement d'un 
fluide entre deux transducteurs disposes en des points 
de mesure espaces dans la direction d'un ecoulement 
du fluide selon lequel la valeur de la vitesse d'ecoule- 

20 ment a mesurer est obtenue en combinant une mesure 
des temps de propagation respectifs de chacun des 
deux signaux acoustiques emis entre les deux points 
dans des sens opposes par les transducteurs avec une 
mesure des dephasages acoustiques respectivement 

25 jnduits dans chaque signal acoustique par la propaga- 
tion de chacun desdits signaux acoustiques dans recou- 
pment, la mesure des temps de propagation respectifs 
de chacun des deux signaux acoustiques emis entre les 
deux points de mesure etant effectuee conformement 

30 au procede selon la presente invention. 

[0025] D'autres caracteristiques et avantages de In- 
vention ressortiront de la description suivante de modes 
particuliers de realisation, faite en reference aux des- 
sins annexes, sur lesquels : 

35 

la Figure 1 est une representation, sous la forme 
d'un schema-bloc, de I'ensemble d'un dispositif de 
comptage de gaz connu auquel est applicable la 
presente invention ; 
^0 - la Figure 2 est un diagramme montrant le principe 
d'une methode connue de determination du temps 
de propagation d'un signal acoustique dans un mi- 
lieu fluide ; 

la Figure 3 est un diagramme agrandi d'une partie 
45 de la Figure 2 montrant le principe d'une methode 
connue de determination du temps de propagation 
d'un signal acoustique par comparaison du signal 
regu avec une valeur de seuil ; 
la Figure 3a est un diagramme analogue a celui de 
50 la Figure 3 montrant la fagon dont on opere sur un 
signal echantillonne, selon une methode connue de 
determination du temps de propagation d'un signal 
acoustique par passage a zero ; 
la Figure 4 est un ensemble de deux diagrammes 
55 analogues a celui de la Figure 3 montrant les ris- 
ques d'erreur dans une methode connue de deter- 
mination du temps de propagation d'un signal 
acoustique faisant appel a une comparaison des si- 
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gnaux regus avec une valeur de seuil ; 
la Figure 5 est un diagramme montrant le principe 
de mise en oeuvre du procede selon I'invention pour 
la determination du temps de propagation d'un si- 
gnal acoustique dans un milieu fluide ; et 5 
la Figure 6 est un diagramme comportant des fa- 
milies de courbes permettant de choisir des valeurs 
de parametres utiles pour la mise en oeuvre du pro- 
cede selon I'invention. 

la Figure 7 est un organigramme montrant un exem- 10 ■■ 
pie des differentes etapes du procede selon inven- 
tion. 



[0026] Differents types de circuits electroniques peu- 
vent etre utilises pour assurer un comptage de fluide par 
une methode de mesure a ultrasons de la vitesse 
d'ecoulement du fluide. La Figure 1 montre le schema 
de~principe d'un exemple de tels circuits electroniques 
qui sont associes a deux transducteurs ultrasonores 1, 

2 disposes a distance Tun de I'autre dans une conduite 

3 dans laquelle s'ecoule un fluide tel qu'un gaz. Les deux 
transducteurs 1 , 2 sont relies a un bloc de commutation 

4 qui comprend deux commutateurs 5, 6 et permet I'uti- 
lisation de chaque transducteur alternativement comme 
emetteur et recepteur. Un module d'emission 14 et un 
module de reception 17 sont relies respectivement aux 
commutateurs 5, 6 du bloc de commutation 4. Le mo- 
dule d'emission 14 comprend un amplificateur opera- 
tionnel 16 et un convertisseur numerique-analogique 
15. Le module de reception 17 comprend au moins un 
amplificateur 18 et un convertisseur analogique-nume- 
rique 19 qui numerise et echantillonne simultanement 
le signal regu. Une source d'energie electrique 7 et un 
module 8 de gestion de I'alimentation electrique sont re- 
lies notamment aux modules d'emission 1 4 et de recep- 
tion 17 ainsi qu'au bloc de commutation 4 et a un micro- 
controleur 10. Le microcontroleur 10 comprend notam- 
ment une horloge a quartz 9, une unite arithmetique et 
logique, des circuits de memoire vive et de memoire 
morte, et peut cooperer avec des circuits d'affichage 1 3, 
une memoire morte reinscriptible 12 et une liaison serie 
11 du type RS232. 

[0027] Les compteurs de gaz destines a equiper cha- 
que abonne d'un reseau de distribution doivent etre a la 
fois precis, fiables et le meilleur marche possible. Ces 
contraintes imposent d'eviter I'utilisation de composants 
couteux, de mettre en oeuvre de petites piles d'alimen- 
tation de longue duree mais de capacite moyenne pour 
des raisons de securite et d'effectuer des calculs selon 
un procede de mesure assurant a la fois une precision 
et une fiabilite elevees tout en restant suffisamment sim- 
ple pour etre econome en energie. 
[0028] Le procede de mesure selon I'invention qui se- 
ra explicite plus loin est ainsi avantageusement appli- 
que a un dispositif de mesure de la vitesse d'ecoulement 
d'un fluide gazeux, permettant d'assurer le comptage de 
la consommation de ce fluide, en mettant en oeuvre 
deux transducteurs ultrasonores disposes a distance 
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I'un de I'autre a Tinterieur d'une conduite, dans le sens 
d'ecoulement du fluide. A titre d'exemple, les transduc- 
teurs ultrasonores peuvent fonctionner a une frequence 
acoustique Fa de I'ordre de 40 kHz et la frequence 
d'echantillonnage Fe des signaux acoustiques regus est 
avantageusement de 320 kHz, c'est-a-dire huit fois plus 
elevee que la frequence Fa. 

[0029] La frequence d'echantillonnage Fe est avanta- 
geusement egale a un multiple de la frequence acous- 
tique Fa. Toutefois, la frequence Fe peut egalement 
prendre d'autres valeurs. Lorsque les valeurs (multiples 
de Fa ou non multiples) de Fe sont trop basses, il est 
preferable d'utiliser par exemple un algorithme d'ap- 
proximation au voisinage du maximum (en valeur abso- 
15 lue) de chaque periode du signal acoustique, voire d'ef- 
fectuer plusieurs mesures successives en dephasant a 
chaque nouvelle mesure le signal d'echantillonnage par 
rapport au signal acoustique considere, ceci afin d'ob- 
tenir une plus grande precision sur les valeurs des am- 
20 plitudes maximales (en valeur absolue). Ces methodes 
seront detaillees ulterieurement. 
[0030] Lorsque les valeurs de Fe sont suffisamment 
elevees I'emploi de telles methodes n'est pas necessai- 
re. 

25 [0031] Le signal acoustique emis par chaque trans- 
ducteur est constitue par exemple par une impulsion 
emise a la frequence acoustique determinee Fa. Cette 
impulsion est par exemple rectangulaire. 
Plusieurs impulsions pourraient egalement etre emises. 

30 [0032] D'une maniere generale, le debit de fluide 
dans une conduite peut s'ecrire : 

cD = (SL72) (Tu - Td) / Tu.Td 

35 

OU 

S designe la section moyenne de passage offerte a 
I'ecoulement entre les deux transducteurs acous- 
40 tiques, 

L designe la distance separant les transducteurs, 
Td est le temps de propagation du signal acoustique 
emis dans le sens de I'ecoulement par le premier 
transducteur (amont) jusqu'a la reception de ce 
45 signal acoustique par le second transducteur 

(aval), 

Tu est le temps de propagation du signal acoustique 
emis dans le sens contraire de I'ecoulement par 
le second transducteur, jusqu'a la reception de ce 
50 signal acoustique par le premier transducteur 

(amont). 

[0033] Si la mesure de ce debit fait en outre intervenir 
une mesure des dephasages acoustiques induits dans 
55 chaque signal acoustique par la propagation de chacun 
des signaux acoustiques emis par les transducteurs 
dans I'ecoulement, le debit de fluide dans une conduite 
peut s'ecrire : 



